Vacuum-polarization corrections to the parity-nonconserving 
6s — 7s transition amplitude in 133 Cs. 

W. R. JohnsonB, I. Bednyakov and G. Soff 
Institut fur Theoretische Physik, Technische Universitdt Dresden, Mommenstrasse 13, D-01062 

Dresden, Germany 



Abstract 

The dominant one- loop radiative corrections to atomic wave functions, those 
associated with vacuum polarization in the nuclear Coulomb field, are eval- 
uated for the 6s — 7s parity-nonconserving (PNC) transition amplitude in 
133 Cs. These corrections increase the size of the PNC amplitude by 0.4% 
and, correspondingly, increase the difference between the experimental value 
of the weak charge Qvk( 133 Cs) and the value predicted by the standard model. 

PACS numbers: 32.80.Ys, ll.30.Er, 31.30.Jv 
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I. INTRODUCTION 



Parity non- conservation (PNC) in atoms, described in the standard model of the elec- 
troweak interaction by exchange of Z bosons between bound electrons and nuclear quarks, 
leads to nonvanishing electric-dipole matrix elements between atomic states with the same 
parity. The nuclear spin-independent part of PNC matrix elements (arising from the vector 
nucleon current) is proportional to a conserved weak charge Qy/i which is sensitive to physics 
beyond the standard model. 

Measurements of the 6s — 7s PNC amplitude in 133 Cs, following the procedure described 
by Bouchiat and Bouchiat [[![, were carried out to 2% accuracy by Gilbert and Wieman ]2[ 
and, more recently, to 0.3% accuracy by Wood et al. @. Following the recent measurements, 
there was a revival of interest in the associated atomic structure calculations of the PNC 
amplitude Indeed, Bennett and Wieman || analyzed differences between experimental 
and theoretical values of amplitudes for allowed transitions, polarizabilities, and hyperfme 
constants in Cs and concluded that the error in the atomic structure calculations of the 
PNC amplitude should be reduced from 1%, the value given in Refs. [§J^]], to 0.4%. Using 
this revised estimate of the accuracy of the calculations, they showed that the experimental 
value of the weak charge of 133 Cs differed from the standard-model value 0, 

Q w ( 133 Cs) = -73.09 ± 0.03, (1) 

by 2.3 a. This difference between experiment and theory, being one of the two largest 
differences reported in the current review of particle physics [J7|, suggested the existence of 
a second neutral Z' boson |8|-|T0"f. Implications of this difference for new physics were also 



reviewed in |TTJ and discussed in |12|Jil 

Breit corrections to the PNC amplitude, which were ignored in and underestimated 
in [§], were shown to decrease the size of the calculated PNC amplitude by a 0.6% by 



Derevianko ||14|| ; this finding was confirmed in Refs. [[y],|T(J. Including Breit corrections 
reduces the difference with the standard model to 1.6 a if the 0.4% error in the calculations 
determined in M is assumed or to 0.9 a if the more conservative 1% error given in Refs. p],|5[] 
is assumed. 

Radiative corrections to PNC matrix elements (w\Hpnc\v) i n the strong Coulomb field of 
a highly-charged one-electron ion were considered recently by Bednyakov et al. |17]]. These 
corrections were decomposed into two parts, radiative corrections to the operator H PNC , and 
radiative corrections to wave functions \v) and \w). The former corrections were evaluated 
for 133 Cs in Refs. [fT8|-21| and are already included in the theoretical value Qw ( 133 Cs) given 
in Eq. ([!]). The dominant part of the wave function radiative correction, which arises from 
the one-loop diagrams Fig. [I], was evaluated for 2s — 2pi/ 2 matrix elements in highly-charged 



one-electron ions in Ref. |17| . Sushkov, in his analysis of the Breit corrections p2"| |, suggested 
that the residual radiative corrections to PNC amplitudes in many-electron atoms could be 
of the same order of magnitude as the Breit corrections. We examine this suggestion further 
in the present paper and find that the one-loop wave-function radiative corrections increase 
the size of the PNC amplitude in 133 Cs by 0.4% and, correspondingly, increases the difference 
between the theoretical and experimental weak charge. Moreover, we find that the one-loop 
wave function correction is insensitive to electron-electron correlation effects. 
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II. CALCULATION 



One-loop radiative corrections (vacuum polarization corrections) in an external Coulomb 
field were considered by Wichmann and Kroll ||23|| , who showed that these corrections simply 
modify the electron-nucleus Coulomb interaction at short range. The modification of the 
Coulomb interaction is described, to leading order in powers of Za, by the Uehling potential: 



El 



We use atomic units here; a = 1/137.036 ... is the fine structure constant and c = a' 1 is 
the speed of light. 

In applications to atomic PNC, the Coulomb potential is also modified at short range 
by finite nuclear size effects, which are described by a Fermi-type charge distribution p(r): 

P ^ l + exp[(r-c nuc )/a nuc ]' ^ 

For 133 Cs, we assume that the central radius c nuc = 5.6748 fm and that the 10%-90% falloff 
distance is t nuc = 2.3 fm, corresponding to a nuc = 0.523 fm and to a root-mean-square radius 



of the nuclear charge distribution -R rms = 4.807 fm. As shown by Fullerton and Rinker |25 
the Uehling potential can be generalized as follows to accommodate the charge distribution 
p{r) of an extended nucleus: 



SV(r) = / dxxp(x) / dtVt 2 - 1 x 
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Corrections to Dirac-Coulomb energies for Is, 2s and 2p electrons in one-electron ions 
obtained by solving the Dirac equation in the composite nuclear + Uehling potential dis- 
cussed above are in found to be in close agreement with sums of the Uehling, Uehling-finite 



nuclear size, and finite nuclear size corrections obtained perturbatively in ||26|| . A comparison 



of the present results with those from |26j are given in Table |. Uehling corrections to ns 



and npi/2 levels of Cs (the levels of primary interest in PNC calculations in Cs) are about 
twice as large as finite nuclear size corrections and have opposite signs as shown in the table 
for levels with n—1 and 2. 

We carry out two calculations of the 6s — 7s PNC amplitude in the modified potential, 
both leading to precisely the same relative correction to the PNC amplitude. The first of 
these calculations, is done at the "weak" Dirac-Hartree-Fock (DHF) level of approxima- 
tion. The perturbation S(j)^ F to a valence electron wave function induced by the weak 
interaction /ipnc satisfies the inhomogeneous DHF equation 

(h + V™ - e f ) 50f = -/ipncC- (5) 

In this equation, Vhf is the HF potential of the closed xenon-like core and e^ F is the eigen- 
value of the unperturbed DHF equation. The perturbed DHF equations are solved to give 
the #0g s F and SffiJ . The PNC amplitude is then given by the sum of two terms: 
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where D is the dipole operator. This calculation leads to a value of the PNC amplitude 
that is 20% lower than the final correlated value found in Refs. p|,|5|,p!6| . In the top panel of 
Table [TI], we show each of the two terms making up the sum in Eq. (Q) determined with and 
without the modified Uehling potential 5V from Eq. @. The one- loop correction is seen to 
increase the size of each terms and their sum by 0.41%. 

The second calculation is done at the "weak" random-phase approximation (RPA) level 
of approximation, in which the class of correlation corrections associated with weak pertur- 
bations of the core orbitals are included in the calculation. This calculation leads to a value 
of the PNC amplitude that is 3% larger than the final correlated value of the amplitude 
given in [f|,[5|,|16| . This class of correlation corrections is included by solving 



(ho + V« F - 6f 



^PNC + VpNC 



/>HF 



(7) 



where V^f c is the weak correction to the HF potential. The resulting PNC amplitude is 
given by 
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In the lower panel of Table |TJ, we give values of the two terms making up the sum in Eq. (^|) 
with and without the Uehling corrections. Again, the one-loop correction is seen to increase 
the size of each of the two terms and their sum by 0.41%. 

From the two calculations above, it is apparent that the short range vacuum-polarization 
corrections are independent of electron-electron correlation. The situation is similar to that 
found for the short range nuclear "skin" correction arising from the difference between 
the neutron and proton radius of the nucleus |]27|j28|] , which also leads to a correlation- 
independent correction to the PNC amplitude. 



III. CONCLUSION 

One-loop radiative corrections to the electron wave functions from vacuum-polarization 
in the nuclear field are evaluated. The resulting wave functions are used to calculate the 
PNC amplitude for the 6s — 7s PNC transition in 133 Cs, leading to a 0.4% increase the size 
of the amplitude. This correction is found to be correlation independent. 

An average of the three most accurate calculations [3fl|l6l , of the 6s — 7s amplitude, 
taking account of Breit corrections and one-loop radiative corrections gives 

E PNC = -9.057 ±0.037 iea x l()- 12 (-Q w /N), 

were we use the estimate from |J for the error in the calculations. Combining this calculated 
amplitude with the experimental value of the amplitude from || leads to an experimental 
weak charge 

Q e vT( 133 Cs) = -72.12 ± (0.28) cxpt ± (0.34) theor 
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This value differs by 2.2 a from the standard model. If we assume a 1% error in the 
theoretical amplitude, then the theoretical component of the error in Q e w vt is increased to 
±(0.74) theor and the difference with the standard model becomes 1.2 a. 

We have considered only the dominant one-loop wave-function radiative corrections in 
the above calculation and ignored the higher-order aZ corrections to the Uehling potential 



discussed by Wichmann and Kroll p3| . The Wichmann-Kroll corrections can be easily 
estimated and are expected to have a negligible effect on the present result. Moreover, the 
"vertex" radiative correction to the PNC amplitude, which contributes about 0.1% to the 
Standard Model value of Qw, is evaluated in Refs. [0,0 using free-particle propagators. 
This radiative correction should be redone using Coulomb-field propagators. If Coulomb- 
field effects were to change the vertex correction by 100% (probably a gross overestimate) 
the Standard Model value of Qw would change by only 0.1%. We therefore expect that 
further changes to the value of Qw from the strong-field corrections to the operator i^PNc 
discussed in [|T7| will be smaller than 0.1%. 
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FIGURES 
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FIG. 1. One loop wave function corrections. The double line represents the electron in the field 
of the nucleus, the wavy lines represent Z bosons or photons and the open triangle represents the 
Z-nucleus vertex. 
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TABLES 



TABLE I. Differences between the present one-electron Dirac energies for Cs (Z=55) in the 
finite nucleus plus Uehling potential and Dirac-Coulomb energies (SE) are compared with pertur- 
bative values of Uehling (Uehl), Uehling-finite nuclear size (Uehl-FS), and finite nuclear size (FS) 
corrections to one-electron energy levels given in J26|. Units: amc 2 (aZ) 4: /(7m 3 ) 



Present 






Ref. || 




State 5E 


Uehl 


Uehl-FS 


FS 


Total 


lsi/2 -0.1419 


-0.2584 


0.0010(1) 


0.1159 


-0.1415(1) 


2s 1/2 -0.1554 


-0.2901 


0.0011(1) 


0.1339(1) 


-0.1551(2) 


2 Pl/2 -0.0100 


-0.0145 


0.0000(1) 


0.0045 


-0.0100(1) 


2p 3 /2 -0.0016 


-0.0016 


0.0000 


0.0000 


-0.0016 


TABLE II. One-loop corrections to the weak HF-PNC amplitude for the 6s- 


— is ri>i Li ampntude 


for 133 Cs are shown in the upper panel and to the weak RPA-PNC amplitude 


in the lower panels. 


Values are listed without and with corrections from Eq. (|j 


). Units: ieao x 10 


~ 12 (-Qw/N) 


Type 


(<p 7s \D\5(f) Gs ) 


( 


5<f) 7s \D\(f) 6s ) 


Epnc 


Weak HF approximation: 










DHF 


2.749183 




-1.014387 


-7.394685 


DHF+ 5V 


2.760414 




-1.018581 


-7.425391 


A(%) 


0.41 




0.41 


0.41 


Weak RPA approximation: 










RPA 


3.457036 




-1.272562 


-9.268581 


RPA+ SV 


3.471169 




-1.277834 


-9.307166 


A(%) 


0.41 




0.41 


0.41 
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